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ABSTRACT 

Context. The nearby Chamaeleon molecular cloud complex is a good laboratory to study the process of low-mass star formation 
since it consists of three clouds with very different properties. Chamaeleon III does not show any sign of star formation, while star 
formation has been very active in Chamaeleon I and may already be finishing. 

Aims. Our goal is to determine whether star formation can proceed in Cha III by searching for prestellar cores, and to compare the 
results to our recent survey of Cha I. 

Methods. We used the Large APEX Bolometer Array (LABOCA) to map Cha III in dust continuum emission at 870 /jm. The map 
is compared with a 2MASS extinction map and decomposed with a multiresolution algorithm. The extracted sources are analysed by 
carefully taking into account the spatial filtering inherent in the data reduction process. 

Results. 29 sources are extracted from the 870 yum map, all of them being starless. The estimated 90% completeness limit is 0.18 M G . 
The starless cores are found down to a visual extinction of 1.9 mag, in marked contrast with other molecular clouds, including Cha I. 
Apart from this difference, the Cha III starless cores share very similar properties with those found in Cha I. They are less dense than 
those detected in continuum emission in other clouds by a factor of a few. At most two sources (< 7%) have a mass larger than the 
critical Bonnor-Ebert mass, which suggests that the fraction of prestellar cores is very low, even lower than in Cha I (< 17%). Only 
the most massive sources are candidate prestellar cores, in agreement with the correlation found earlier in the Pipe nebula. The mass 
distribution of the 85 starless cores of Cha I and III that are not candidate prestellar cores is consistent with a single power law down to 
the 90% completeness limit, with an exponent close to the Salpeter value. A fraction of the starless cores detected with LABOCA in 
Cha I and III may still grow in mass and become gravitationally unstable. Based on predictions of numerical simulations of turbulent 
molecular clouds, we estimate that at most 50% and 20% of the starless cores of Cha I and III, respectively, may form stars. 
Conclusions. The LABOCA survey reveals that Cha III, and Cha I to some extent too, is a prime target to study the formation of 
prestellar cores, and thus the onset of star formation. Getting observational constraints on the duration of the core-building phase prior 
to gravitational collapse will be necessary to make further progress. 

Key words. Stars: formation - ISM: individual objects: Chamaeleon III - ISM: structure - ISM: evolution - ISM: dust, extinction - 
Stars: protostars 

in 

1 . Introduction 11998b see Appendix Bl of Belloche et al. 201 1 for more details) 

^— ' . and mainly consists of three molecular clouds, Chamaeleon I, 

, W ith the advent of large (sub)mm bolometer arrays, the search n and m (he reafter Cha I, II, and III), that have very different 
for cold dense cores in molecular clouds is becoming more ef- degr ees of star formation activity. Their populations of young 
ficient. These dense cores correspond to the earliest stages of stellar ob j ects (Y SOs) have been well studied from the infrared 
the birth of stars and their study is essential to understand the to x . rays; nearly one order of magn i tu de more YSOs have been 
process of star formation. In particular, unbiased searches for found in cha T comp ared to Cha II - while Cha III do es not seem 
prestellar coresU and protostars are needed to better understand to contain any Y SO observable at these wavelengths (iPersi et all 
their formation and derive their lifetimes. |2003l: iLuhmad 120081) . The three clouds have similar masses as 

In this respect, the Chamaeleon molecular cloud complex traced with 13 CO 1-0 (~ 1000 M Q ) but the fraction of denser 
is of particular interest. It is one of the nearest low-mass star gas traced with C ls O 1-0 is the hi ghest in Cha I (24%) whil e 
forming regions (150-180 pc. lWhittet et al.ll997t lKn ude & Hflgl it i s the lowest in Cha III (4%) (iMizuno etalJll999l l200lh . 

Finally, several indications of jets/outflows were found in Cha I 

* Based on observations carried out with the Atacama Pathfinder (Mattila et al. 1989: Gomez et alJl2004l:IWang & Henning||2006l : 
Experiment telescope (APEX). APEX is a collaboration between iBelloche e t alj2006). Only three Herbig-Haro objects are know n 
the Max-Planck Institut fur Radioastronomie, the European Southern in Cha II and none has been found in Cha III dSchwartzlll977l) . 
Observatory, and the Onsala Space Observatory. Therefore, Cha III is the least active region in the Chamaeleon 

** The FITS file of Fig. [2] is available in electronic form at the CDS cloud complex, 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5). 

1 A prestellar core is usually defined as a starless core tha t is gravi- Among the three clouds of the Chamaeleon complex, Cha III 

tationally bound dAndre et al.ll200Cl ; ldi Francesco et alj|2007h . Starless stands out with a prominent filamentary structure as revealed 
means that does not contain any young stellar object. by the InfraRed Astronomical Satellite (IRAS) at 100 //m (see 



- 1—1 

X 



1 



A. Belloche et al.: Will the starless cores in Chamaeleon I and III turn prestellar?' 



Figs. 5 and 7 of iBoulanger et al.l fl998). It consists of a system 
of interwoven filaments that can be disentangled wit h velocity 
inform ation derived from molecular line emission (Gah m et alj 
2002). With an angular resolution of 45", better than that of 
the 13 CO/C ls O surveys mentioned above, 38 small clumps em- 
bedded in these filaments were detected in C 18 1-0 with the 
Swedish ESO Submillimeter Telescope (SEST) in the course of 
a survey targetting the col umn density peaks of cold dust emis- 
sion based on IRAS data dGahm et alj|2002l) . The clumps have 
mean densities in the range 1-8 xlO 4 cirT 3 and their internal ki- 
netic energy is dominated by turbulence. Most of them are far 
from virial equilibrium, suggesting that they are not sites of cur- 
rent star formation. 

Since Cha III was mapped in the 1-0 transition of CO and its 
isotopologues, a transition tracing low densities and a molecule 
suffering from depletion at high density, only the distribution of 
its low d ensity gas is relatively well known. However, in contrast 
to Cha I ( Bello che et al.ll20lH hereafter Pap er I) and in a limite d 
way Cha II (low-sensitivity 1.3 mm survey of lYoung et al.l2 005). 
no systematic survey for (sub)mm dust continuum emission has 
been performed in Cha III. Therefore, the population of dense, 
prestellar cores in Cha III is completely unknown. On the one 
hand, the absence of signposts of active star formation in the 
protostellar and pre-main-sequence phases in Cha III could be 
the result of environmental conditions different from those pre- 
vailing in, e.g., Cha I, and unfavorable to star formation. On the 
other hand, star formation could just be starting in Cha III and 
finding prestellar cores would favor this interpretation. 

To unveil the present status of the earliest stages of star for- 
mation in Cha III and compare this cloud to Cha I, we car- 
ried out a deep, unbiased dust continuum survey of Cha III at 
870 pm with the Large APEX Bolometer Camera at the Atacama 
Pathfinder Experiment (APEX). The observations and data re- 
duction are described in Sect. [2] Section[3]presents the maps and 
the method used to extract the detected sources. The properties 
of these sources are analysed in Sect. [4] The implications are 
discussed in Sect. [5] Section [6] gives a summary of our results 
and conclusions. 



2. Observations and data reduction 

2.1. Extinction map from 2MASS 

We derived an extinction map of Cha III from the publicly avail- 
able 2MASS0 point source catalog in the same way as we did 
for Cha I (see Paper I), except that no source filtering was ap- 
plied since there are no known embedded YSOs in Cha III. We 
used the same resolution of 3' (FWHM) with a pixel size of 1 .5'. 
With these parameters, most pixels have at least 10 stars within 
a radius equal to FWHM/2. Only a few pixels contain fewer 
stars, the minimum being 5 stars for 6 pixels. The resulting map 
is shown in Fig.Q] The typical rms noise level in the outer parts 
of the map is 0.4 mag, corresponding to a 3cr detection level of 
1 .2 mag for an FWHM of 3'. This rms noise level is however ex- 
pected to increase towards the higher-extinction regions because 
of the decreasing number of stars per element of resolution. 



2 The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 
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Fig. 1. Extinction map of Cha III derived from 2MASS in radio 
projection. The projection center is at (a, #)j2000 = (12 h 42 m 24 s , 
-79°43'48"). The contours start at Ay = 3 mag and increase 
by steps of 1.5 mag. The dotted lines are lines of constant right 
ascension. The angular resolution of the map (HPBW = 3') is 
shown in the upper right corner. The five fields selected for map- 
ping with LABOCA are delimited with dashed lines. The field 
of view of LABOCA is displayed in the lower right corner. 



2.2. 870 pm continuum observations with APEX 

The region of Cha III with a visual extinction higher than 3 mag 
was selected on the basis of the extinction map derived from 
2MASS as described above (Sect. 12. U . It was divided into five 
contiguous fields labeled Cha3-North, East, Center, South, and 
West with a total angular area of 0.93 deg 2 (see Fig . |T} . The five 
fields were mapped in continuum e mission with the La rge APEX 
BOlometer CAmera (LABOCA, Sirin go etaT] 12009) operating 
with about 250 working pixels in the 870 pm atm ospheric win- 
dow at the APEX 12 m submillimeter telescope (Gtiste n et al.l 
120061) . The central frequency of LABOCA is 345 GHz and its 
angular resolution is 19.2" (HPBW). The observations were car- 
ried out for a total of 88 hours in September and November 2010, 
under excellent (t 870 '™ =0.12) to moderate (t 87 ^ = 0.43) at- 

v zenith y v zemtn ' 

mospheric conditions. The sky opacity was measured every 1 
to 2 hours with skydips. The focus was optimised on rj Carina, 
Mars, or G34.26+0.15 at least once per day/night. The pointing 
of the telescope was checked every 1 to 1 .5 hour on the nearby 
quasar PKS 1057-79 and was found to be accurate within 2.5" 
(rms). The calibration was performed with the secondary cali- 
brators IRAS 13134-6264, G5.89-0.39, G34.26+0.15, or NGC 
20 71 that were obse rved every 1 to 2 hours (see Table A.l 
of Siringo et al. 2009). Measurements on the primary calibrator 
Mars were also used. 

The observations were performed on the fly with a rectangu- 
lar pattern ("OTF"). The OTF maps were performed with a scan- 
ning speed of 2 arcmin s _1 and were alternately scanned in right 
ascension and declination, with a random position angle between 
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-12° and +12° to improve the sampling and reduce striping ef- 
fects. 

2.3. LABOCA data reduction 

The LABOCA data were reduced with the BoA softwar^l fol- 
lowing the iterative procedure described in Paper I. The only dif- 
ference is that the baselines of the individual OTF maps scanned 
in declination had to be removed subscan-wise rather than scan- 
wise to improve the flatness of the background level. The grid- 
ding was done with a cell size of 6. 1" and the map was smoothed 
with a Gaussian kernel of size 9" (FWHM). The angular reso- 
lution of the final map is 21.2" (HPBW) and the rms noise level 
is 11.5 mJy/21.2"-beam (see Sect. EDI 

The spatial filtering properties of the data reduction and the 
convergence of the iterative process are analysed in AppendixlAl 
in the same way as was done for Cha I (Paper I). In short, the 
Cha III dataset that was obtained with rectangular scanning pat- 
terns seems to suffer slightly more from spatial filtering than the 
Cha I dataset that combined rectangular and more compact spi- 
ral scanning patterns. The reasons for this counterintuitive result 
are unclear. It may be related to the slightly smaller number of 
well- working pixels for the Cha III dataset compared to the Cha I 
one. 



3. Basic results and source extraction 

In the following, we make the same assumptions as in Paper I 
(see its Appendix B) to derive the physical properties of the 
detected sources, in par ticular we assu me a uniform dust tem- 
perature of 12 K (s ee iToth et al.l l2000h . a distance of 150 pc 
(Whittet et al. 1997), a dust mass opacity, k%^, of 0.01 cm 2 per 
gram of gas+dust, a gas-to-dust mass ratio of 100, and a mean 
molecular weight per free particle, fi, of 2.37. These assumptions 
are not repeated in the following, except in the few cases where 
there could be an ambiguity. 



3.1. Maps of dust continuum emission in Cha III 

The final 870 fim continuum emission map of Cha III obtained 
with LABOCA is shown in Fig. [2] Pixels with a weight (1/cr 2 ) 
smaller than 3500 beam 2 /Jy 2 are masked. The resulting map 
contains 0.32 megapixels (out of 0.76 that contain some signal), 
corresponding to a total area of 0.92 deg 2 (6.3 pc 2 ). The mean 
and median weights are 5717 and 5798 beam 2 /Jy 2 , respectively. 
The noise distribution is fairly uniform and Gaussian. The av- 
erage noise level is 11.5 mJy/21.2"-beam. This translates into 
an H2 column density of 1.0 x 10 21 cm 2 , and corresponds to 
a visual extinction Ay ~ 1.1 mag with Ry = 3.1 (see the other 
assumptions in Appendix B of Paper I). 

The dust continuum emission map of Cha III reveals many 
weak, spatially resolved sources. In contrast to Cha I, not a single 
unresolved, compact source is detected, which is consistent with 
the absence of signposts of star formation at other wavelengths 
in Cha III. Figure [3] presents all the detected structures in more 
detail. The most prominent one is the dense core in the north- 
ern part of field Cha3-North (Fig. [3^). All other detected struc- 
tures are much fainter. Although there is no direct detection of 
any large-scale, filamentary structure (see Fig. HJ}), the distribu- 
tion of detected sources is much reminiscent of the filamentary 



structures seen in the far-infrared (see Sect. [TJ. In field Cha3- 
East, most detected sources (Cha3-C21, 7, 17, 12, and fainter 
3<x compact structures) are distributed along a 1.7 pc long fila- 
ment (Fig.[3l>). A second, shorter (0.3 pc), filamentary structure 
that may connect to the latter is suggested by the spatial distri- 
bution of Cha3-C19, 8, and a fainter 3cr compact structure to 
the south-west. In field Cha3-Center, Cha3-C3, 4, 13, 5, 16, and 
22 are remarkably aligned and nearly equally distributed along 
a straight line of length 0.9 pc (Fig. [3]:). Finally, the sources de- 
tected in field Cha3-South also suggest the existence of a fila- 
ment of length 0.8 pc (Fig.[3]d). 

Even if the filamentary structure of Cha III is not directly 
seen with LABOCA, we expect that it will be detected with 
the Her schel Space Obse rvatory in the frame of the Gould Belt 
Survey (I Andre et al.ll2010l) . Filaments have been detected in all 
clouds analysed from this survey so far, and a close connection 
between these fila ments and the formation of dense cores has 
been established (lAndre etal.ll20TotlMen , shchikovetalJl2010l: 
lArzoumanian et al.ll201 ll) . 

3.2. Masses traced with LABOCA and 2MASS 

The total 870 fim flux in the whole map of Cha III is about 
42.7 Jy. This translates into a cloud mass of 22.6 M . It cor- 
responds t o 1.2% of the tota l mass traced by CO in Cha III 
( 1 890 IVL. lMizuno et al.l200lt) 2.0% of the mass traced by 13 CO 
(1 100 Mm. iMizuno et alJll999h. a nd 54% of the mass traced by 
C ls O (42 IVU lMizuno et al]1999h . In Cha I, these fractions were 
5.9%, 7.7%, and 27-32%, respectively (Paper I). 

The extinction map shown in Fig.Q~]traces larger scales than 
the 870 /mi dust emission map. The median and mean extinc- 
tions over the 0.92 deg 2 covered with LABOCA are 2.8 and 
3.0 mag, respectively. Although our survey was designed to 
cover the regions above 3 mag, a significant fraction of the re- 
sulting map that is based on rectangular fields comprises regions 
below 3 mag, which explains these low median and mean extinc- 
tions. Assuming an extinction to H2 column density conversion 
factor of 9.4 x 10 20 cirT 2 mag -1 (for Ry = 3.1, see Appendix B. 3 
of Paper I), we derive a total gas+dust mass of 401 M . However, 
96% of this mass, i.e. 386 M , is at Ay < 6 mag. With the ap- 
propriate conversion factor for Ay > 6 mag (see Appendix B 
of Paper I), the remaining mass is reduced to 9.3 M Q , yield- 
ing a more accurate estimate of 395 M for the total mass of 
Cha III traced with the extinction. It is much lower than the 
masses traced by CO and 13 CO mentioned above. Since the 
latter mass values (at least those derived from CO emission) 
were based on integrations over much larger areas where CO 
and 13 CO still emi t significantljQ (see region VI in Fig. 2 of 
iMizuno et al .1 120011) . we consider the mass derived from the ex- 
tinction map as the best estimate to compare with. Thus the mass 
traced with LABOCA represents about 5.7% of the cloud mass, 
which is slightly less than in Cha I (7.5%, see Paper I). Given 
that the median extinction is 2.8 mag in the extinction map, i.e. 
about 2.5 times the rms sensitivity achieved with LABOCA, the 
missing 94% were lost not only because of a lack of sensitiv- 
ity but also because of the spatial filtering due to the correlated 
noise removal (see Sect. LA.lt . Finally, we estimate the average 
density of free particles in the field covered with LABOCA. We 
assume that the depth of the cloud along the line of sight is equal 
to the square root of its projected surface, i.e. 2.5 pc. This yields 



3 See http://www.mpifr-bonn.mpg.de/div/submmtech/software/boa/ 
boa_main.html. 



4 This ca veat does not concer n the mass derived from C I8 1-0 since 
the map of Miz uno et alj Q999;) has more or less the same size as the 
LABOCA map. 
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Fig. 3. Detailed 870 /urn continuum emission maps of Cha III ex- 
tracted from the map shown in Fig. [2] The flux density greyscale 
is shown on the right of each panel and labeled in Jy/21.2"- 
beam. It has been optimized to reveal the faint emission with 
a better contrast. The contours start at a and increase with a 
step of a, with a = 34.5 mJy/21.2"-beam, i.e. 3 times the rms 
noise level. The dotted blue contour corresponds to -a. The an- 
gular resolution of the map is shown in the lower left corner of 
each panel (HPBW = 21.2"). The white plus symbols and el- 
lipses show the positions, sizes (FWHM), and orientations of 
the Gaussian sources extracted with Gaussclumps in the filtered 
map shown in Fig. [4^. The sources are labeled like in the first 
column of Table [3] The (red) crosses show the p eak position 
of the 38 clumps detected in C ls O 1-0 with SEST dGahm et alj 
2002). a Field Cha3-North. 



an average density of ~ 420 cm" 3 . This value is very similar to 
the average density estimated for Cha I (380 cm" 3 , see Paper I). 
Alternatively, if we assume that Cha III is filamentary and that 
its depth is rather similar to its typical minor size in the plane of 
the sky (roughly 1 pc), then its average density becomes 1 . 1 x 10 3 
cirT 3 . 

The average density derived above corresponds to an aver- 
age thermal pressure P t h/^B of 0.5-1.3 x 10 4 cm" 3 K, assuming 
a kinetic temperature of 12 K. The C ls O and CO 1-0 linewidths 
measured in Cha III are on the order of 0.9 and 2.6 km s , 
respectively dBoulanger et alj|1998l iMizuno et ail 1 19991 1200 lh . 
The 13 CO 1-0 emission tra ces gas of densities on the order 
of 1 x 10 3 cm 3 in Cha III dMizuno et all 1 19981) and its typi- 
cal linewidth must be between those of C 18 and CO 1-0. We 
take 1 km s" 1 as a representative valu e, which is the typical 
value measured by Miz uno et al.l (1 1998b in the small clouds of 
the Chamaeleon complex. The turbulent pressure is defined as 
^turb = A"WH«cr?. with fi the mean molecular weight per free 
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Fig. 3. (continued) b Field Cha3 -East. 



particle, oth the atomic mass of hydrogen, n the average free- 
particle density, and cr NT the non-thermal rms velocity disper- 
sion derived from the linewidth. We obtain a turbulent pressure 
Aurb/^B °f 2.1-5.6 x 10 4 cm" 3 K. The turbulent pressure domi- 
nates by a factor ~ 4 over the thermal pressure in Cha III. It is on 
the same order as the total pressure of the ISM in the mid-plane 
of the Galaxy (2 x 10 4 cm" 3 K, see Cox 2005). The situation is 
similar in Cha I. 



3.3. Source extraction and classification 

The source extraction from the 870 fim map is performed in the 
same way as for Cha I (Paper I). The map is first decomposed 
into different scales with our multiresolution program based on 
a median filter. The total fluxes measured in the sum maps at 
scales 3 to 7 (see definition in Appendix C of Paper I) are listed 
in Table[T] as well as the corresponding masses. About half of the 
total flux is emitted by structures smaller than ~ 200" (FWHM), 
and only 8% by structures smaller than ~ 60". Table[T]shows that 
the fraction of continuum flux at small scales (< 60") is slightly 
smaller than in Cha I. This may suggest that the structures in 
Cha III are less centrally peaked or it could simply be a bias due 
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Fig. 3. (continued) d Field Cha3-South. 
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Fig. 3. (continued) e Field Cha3-West. 



to the slightly higher sensitivity of the Cha III survey. The sum 
map at scale 5 and its associated smoothed map are shown in 
Fig- El The sum of these two maps is strictly equal to the original 
map of Cha III shown in Fig. [2] 

The sources are extracted with the Gaussian fitting pro- 
gram Gaussclumps with the same parameters as for Cha I (see 
Paper I). The sum maps at scales 1 to 7 were decomposed into 
0, 0, 2, 16, 29, 38, and 39 Gaussian sources, respectively, and 



Table 1. Continuum flux distribution in Cha III and comparison 
to Cha I. 



Scale 


Typical size 


Flux 


Mass 


F/F tot 


F/F tot (Cha I) 






(Jy) 


(M ) 


(%) 


(%) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



3 


< 60" 


3.5 


1.9 


8 


11 


4 


< 120" 


10.3 


5.4 


24 


22 


5 


< 200" 


23.0 


12.2 


54 


49 


6 


< 300" 


36.9 


19.5 


86 


82 


7 


-all 


43.2 


22.8 


101 


99 




all 


42.7 


22.6 


100 


100 



Notes. The last row corresponds to the full map, while rows 1 to 5 
correspond to the sum of the filtered maps up to scale i listed in the first 
column (i.e. the sum map at scale i). Column 2 gives the range of sizes 
of the sources that significantly contribute to the emission with more 
than 40% of their peak flux density (see Col. 2 of Table C. 1 of Paper I). 
Columns 5 and 6 give the fraction of flux detected in each map, for 
Cha III and I, respectively. 



the full map was decomposed into 39 Gaussian sources. These 
counts do not include the sources found too close to the noisier 
map edges (weight < 4400 beam 2 /Jy 2 ), which we consider as 
artefacts. 

We now consider the results obtained with Gaussclumps for 
the sum map at scale 5 (i.e. the map shown in Fig. |4^), which 
is a good scale to characterize sources with FWHM < 120" as 
shown in Appendix C of Paper I. The positions, sizes, orienta- 
tions, and indexes of the 29 extracted Gaussian sources are listed 
in Table [2] in the order in which Gaussclumps found them. We 
looked for associations with sources in the SIMBAD astronomi- 
cal database. We used SIMBAD4 (release 1.171) as of February 
10 th , 2011. None of the 29 Gaussclumps sources is associated 
with a SIMBAD object within its FWHM ellipse. As a result, 
we consider that all these 29 sources are starless cores. 



4. Analysis 

4. 1 . Comparison with the extinction map 

The extinction map derived from 2MASS is overlaid on the 
870 pm dust continuum emission map of Cha III in Fig. [5] The 
overall correspondence is relatively good, most of the contin- 
uum emission coinciding with the peaks of the extinction map. 
The 3cr H2 column density sensitivity limit of the 870 pm map 
is 3.0 x 10 21 crrT 2 , which corresponds to Ay ~ 3.2 mag. Most 
of the continuum emission detected at 870 pm is above the 
contour level Ay — 4.5 mag, and most extended regions with 
3 < Ay < 4.5 mag traced by the extinction map are not detected 
at 870 pm. This is partly due to a lack of sensitivity and to the 
spatial filtering related to the sky noise removal since these low- 
extinction regions have sizes on the order of 5'-10', comparable 
to the field of view of LABOCA. On the other hand, there are 
also a few 870 pm sources detected below 4.5 mag that are not 
seen in the extinction map, most likely because of its poor an- 
gular resolution (e.g. Cha3-C6 and 23 in fields Cha3-Center and 
South, respectively). 

4.2. Comparison with C 18 1-0 

The p eak positions of the 38 clumps detected by Gah m et all 
(2002) in C 18 1-0 emission are shown as (red) crosses in 
Figs. [3^-e. One caveat to keep in mind is that the C 18 1-0 
survey was biased and covered only parts of the field mapped 
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Table 2. Sources extracted with Gaussclumps in the 870 fim continuum sum map of Cha III at scale 5. 





R.A. 


Decl. 


f a 

,/peak 


/lot" 


maj." 


min." 


P.A." 


S" 




(J2000) 


(J2000) 


(Jy/beam) 


(Jy) 


(") 


(") 


O 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


12:54:53.58 


-78:52:24.2 


0.183 


2.928 


116.0 


62.1 


-77.9 


84.8 


2 


12:54:23.02 


-78:51:58.0 


0.116 


1.070 


71.9 


57.7 


55.2 


64.4 


3 


12:44:24.79 


-80:08:50.2 


0.100 


0.380 


56.2 


30.4 


11.2 


41.3 


4 


12:43:35.51 


-80:11:11.7 


0.099 


0.831 


77.5 


48.9 


-20.1 


61.6 


5 


12:42:03.29 


-80:16:06.6 


0.093 


0.699 


67.4 


50.4 


-36.7 


58.3 


6 


12:47:38.53 


-80:07:23.2 


0.089 


0.562 


55.7 


50.9 


21.7 


53.2 


7 


12:52:17.60 


-79:25:42.1 


0.087 


0.434 


59.2 


38.0 


-69.4 


47.4 


8 


12:54:20.65 


-79:30:08.0 


0.085 


0.716 


96.2 


39.4 


52.0 


61.5 


9 


12:53:40.41 


-78:59:34.7 


0.085 


0.607 


86.1 


37.4 


48.5 


56.8 


10 


12:54:08.10 


-78:53:51.0 


0.083 


0.301 


59.8 


27.4 


-76.5 


40.5 


11 


12:52:12.93 


-79:22:45.4 


0.079 


0.376 


56.5 


37.7 


-56.9 


46.1 


12 


12:52:56.46 


-79:39:49.4 


0.076 


0.586 


77.4 


44.7 


-5.1 


58.8 


13 


12:42:48.27 


-80:12:36.8 


0.077 


0.421 


65.3 


37.5 


17.1 


49.5 


14 


12:36:46.21 


-80:34:53.5 


0.075 


0.326 


72.7 


26.8 


-32.8 


44.1 


15 


12:55:53.07 


-79:03:51.3 


0.073 


0.406 


63.0 


39.5 


83.1 


49.9 


16 


12:40:45.58 


-80:19:17.0 


0.069 


0.384 


67.6 


36.9 


-9.7 


50.0 


17 


12:52:49.09 


-79:33:06.9 


0.070 


0.101 


30.4 


21.4 


82.9 


25.5 


18 


12:54:37.43 


-78:53:03.5 


0.071 


0.173 


35.4 


30.7 


-20.6 


33.0 


19 


12:54:40.89 


-79:29:03.3 


0.069 


0.242 


50.0 


31.7 


30.5 


39.8 


20 


12:54:47.60 


-78:53:39.6 


0.070 


0.197 


38.2 


33.2 


54.7 


35.6 


21 


12:51:39.85 


-79:23:29.5 


0.063 


0.211 


66.7 


22.7 


-1.9 


38.9 


22 


12:39:36.53 


-80:22:24.5 


0.065 


0.211 


44.4 


32.9 


-85.1 


38.2 


23 


12:35:44.30 


-80:25:21.6 


0.064 


0.498 


120.0 


29.4 


-23.8 


59.4 


24 


12:36:13.68 


-80:34:09.3 


0.063 


0.432 


67.6 


45.5 


-34.9 


55.5 


25 


12:48:31.42 


-79:38:46.7 


0.061 


0.153 


52.4 


21.4 


-89.9 


33.5 


26 


12:56:07.33 


-79:03:22.2 


0.059 


0.073 


26.2 


21.2 


-79.9 


23.6 


27 


12:38:10.69 


-80:19:35.8 


0.058 


0.135 


49.3 


21.2 


-10.0 


32.3 


28 


12:29:06.05 


-80:10:53.3 


0.058 


0.106 


38.7 


21.2 


-30.6 


28.7 


29 


12:51:36.11 


-79:31:17.1 


0.058 


0.155 


41.2 


29.3 


86.1 


34.7 



Notes. Scale 5 is defined in Sect. l3.3l and TablefJJ (a) Peak flux density (in Jy/2 1.2" -beam), total flux, FWHM along the major and minor axes, and 
position angle (east from north) of the fitted Gaussian. {b) Mean source size, equal to the geometrical mean of the major and minor FWHM. 



with LABOCA. For instance, Cha3-C9, 15, and 26 in field Cha3- 
North, Cha3-C8 and 19 in field Cha3-East, and Cha3-C6 in field 
Cha3-Center were not covered. A second caveat is that the C 18 
maps were undersampled by a factor of ~ 2, with a step of 1' 
in fields Cha3-North and East (but 30" in the brightest regions), 
and 40" in fields Cha3-Center, South, and West. 

About twenty 870 //m sources only are detected in the fields 
covered by the C 18 observations, which implies a detection rate 
with LABOCA a factor of 2 lower and shows that the C ls O 
observations were sensitive to lower density gas, as expected. 
Figures [3^-e show that there is no one-to-one correspondence 
between the peak positions of the 870 fj.m and C ls O sources. 
Only seven C 18 sources peak within 1' from the peak position 
of a 870 jum source. This may partly be due to the undersam- 
pling of the C 18 maps, but most likely results from depletion 
that affects C 18 at high density. This confirms that the 870 fim 
emission traces the high density regions better than the C 18 1-0 
emission. The 870 fim map should therefore give a better census 
of the potential future sites of star formation in Cha III. 

4.3. Properties of the starless cores 

The properties of the 29 starless sources detected with LABOCA 
are listed in Table [3] and their distribution is shown in Figs. [6] 
and [7] The column density (Col. 5) and masses (Cols. 7-9) are 
computed with the fluxes fitted with Gaussclumps or directly 
measured in the sum map at scale 5. As a caveat, we remind 
the reader that the assumption of a uniform temperature may be 



inaccurate and bias the measurements of the masses and column 
densities, as well as the mass concentration (or equivalently the 
density contrast). A dust temperature drop toward the center of 
starless dense cores is possible (see Appendix B.2 of Paper I and 
references therein). 

4.3.1. Extinction 

The visual extinctions listed in Table [3] and plotted in Fig. |6^ 
are extracted from the extinction map derived from 2MASS (see 
Sect. I2.U . Given the lower resolution of this map (HPBW = 
3') compared to the 870 pm map, it provides an estimate of the 
extinction of the environment in which the 870 yt/m sources are 
embedded. 

The 870 fim sources are found down to a visual extinc- 
tion, Ay, of ~ 1.9 mag (as traced with 2MASS at low an- 
gular resolution). This is in marked contrast with the thresh- 
old Ay ~ 5-7 mag above which starless sources are found 
in other lo w-mass star forming r egions lik e Cha I (Paper I) , 
Ophiuchus dJohnstone et al.ll2004l), Perseus dEnoch et al.l l2006: 
iKirk et alJ 120061). Taurus (iGoldsmifh et alJ 120081) . and Aquila 
dAndre et al.ll201 lb . We note however that, with its high sensi- 
tivity, Herschel revealed many starless cores at low extinction 
in the nearby Polaris flar e region, a non-star-forming molecular 
cloud dAndre et al.ll2010t) . 

About half of the Cha III sources are found at Ay < 5 mag 
(median 4.9 mag), and the distribution peaks at Ay ~ 4 mag. 
Both the median and mean as well as the peak of the distribution 
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Table 3. Characteristics of starless sources extracted with Gaussclumps in the 870 pm continuum sum map of Cha III at scale 5. 



Name 


N gd " 


FWHM" 


R* b 


^peak 






M m e 


M 50" e 


Cm 1 




^peak 


n * 

"mean 


"50 


Cn l 






(1000 AU) 2 




(10 21 cm- 2 ) 


(mag) 


(M G ) 


(M ) 


(M Q ) 


(%) 




(10 5 cnr 


3 ) (10 4 cm" 


3 ) (10 4 cm" 3 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


Cha3-Cl 


1 


17.1 x 8.8 


2.0 


16 


5.0 


0.097 


1.55 


0.29 


33( 2) 


1.12 


8.6 


3.0 


19.7 


4.4(0.3) 


Cha3-C2 


2 


10.3 x 8.0 


1.3 


10 


8 


0.061 


0.57 


0.23 


27( 3) 


0.55 


5.5 


2.7 


15.5 


3.5(0.4) 


Cha3-C3 


3 


7.8 x 3.3 


2.4 


8.8 


4.2 


0.053 


0.20 


0.11 


47( 7) 


0.35 


4.7 


5.6 


7.6 


6.2(0.9) 


Cha3-C4 


4 


11.2 x 6.6 


1.7 


8.6 


3.7 


0.052 


0.44 


0.15 


35( 5) 


0.45 


4.6 


2.5 


10.1 


4.6(0.6) 


Cha3-C5 


5 


9.6 x 6.9 


1.4 


8.1 


5.7 


0.049 


0.37 


0.13 


37( 5) 


0.40 


4.3 


2.5 


9.0 


4.9(0.7) 


Cha3-C6 


6 


7.7 x 6.9 


1.1 


7.8 


1.9 


0.047 


0.30 


0.13 


37( 6) 


0.36 


4.2 


2.7 


8.7 


4.8(0.7) 


Cha3-C7 


7 


8.3 x 4.7 


1.8 


7.6 


5.1 


0.046 


0.23 


0.10 


44( 7) 


0.32 


4.1 


3.4 


7.1 


5.8(0.9) 


Cha3-C8 


8 


14.1 x 5.0 


2.8 


7.4 


4.5 


0.045 


0.38 


0.11 


39( 6) 


0.40 


4.0 


2.3 


7.7 


5.2(0.8) 


Cha3-C9 


9 


12.5 x 4.6 


2.7 


7.4 


4.9 


0.045 


0.32 


0.11 


41(7) 


0.37 


4.0 


2.6 


7.3 


5.4(0.9) 


Cha3-C10 


10 


8.4 x 2.6 


3.2 


7.2 


6 


0.044 


0.16 


0.083 


52(10) 


0.30 


3.9 


5.6 


5.7 


6.9(1.3) 


Cha3-Cll 


11 


7.9 x 4.7 


1.7 


7.0 


4.9 


0.042 


0.20 


0.10 


42( 7) 


0.29 


3.7 


3.2 


6.8 


5.5(1.0) 


Cha3-C12 


12 


11.2 x 5.9 


1.9 


6.7 


7 


0.040 


0.31 


0.10 


40( 7) 


0.34 


3.6 


2.1 


6.8 


5.3(1.0) 


Cha3-C13 


13 


9.3 x 4.6 


2.0 


6.8 


4.6 


0.041 


0.22 


0.10 


41(7) 


0.30 


3.6 


2.8 


6.8 


5.4(1.0) 


Cha3-C14 


14 


10.4 x 2.4 


4.3 


6.6 


8 


0.040 


0.17 


0.083 


48( 9) 


0.30 


3.5 


4.8 


5.6 


6.3(1.2) 


Cha3-C15 


15 


8.9 x 5.0 


1.8 


6.4 


5.3 


0.039 


0.21 


0.088 


44( 9) 


0.28 


3.5 


2.6 


6.0 


5.8(1.1) 


Cha3-C16 


16 


9.6x4.5 


2.1 


6.0 


5.4 


0.037 


0.20 


0.082 


44( 9) 


0.27 


3.2 


2.5 


5.6 


5.8(1.2) 


Cha3-C17 


17 


3.3 x 2.1 


1.6 


6.1 


5.7 


0.037 


0.053 


0.048 


77(21) 


0.18 


3.3 


10.7 


3.3 


10.0(2.7) 


Cha3-C18 


18 


4.3 x 3.3 


1.3 


6.2 


9 


0.038 


0.091 


0.12 


30( 5) 


0.21 


3.4 


6.1 


8.4 


4.0(0.7) 


Cha3-C19 


19 


6.8 x 3.5 


1.9 


6.0 


4.8 


0.036 


0.13 


0.091 


40( 8) 


0.23 


3.2 


3.9 


6.2 


5.2(1.1) 


Cha3-C20 


20 


4.8x3.8 


1.2 


6.1 


9 


0.037 


0.10 


0.083 


45( 9) 


0.22 


3.3 


4.8 


5.6 


5.9(1.2) 


Cha3-C21 


21 


9.5 x 2.1 


4.5 


5.5 


5.9 


0.033 


0.11 


0.060 


55(14) 


0.22 


3.0 


4.5 


4.0 


7.3(1.8) 


Cha3-C22 


22 


5.8 x 3.8 


1.5 


5.7 


4.3 


0.034 


0.11 


0.073 


47(10) 


0.21 


3.0 


3.8 


5.0 


6.1(1.4) 


Cha3-C23 


23 


17.7 x 3.1 


5.8 


5.6 


2.8 


0.034 


0.26 


0.069 


49(11) 


0.32 


3.0 


2.4 


4.7 


6.4(1.5) 


Cha3-C24 


24 


9.6 x 6.0 


1.6 


5.5 


7 


0.033 


0.23 


0.091 


37( 8) 


0.27 


3.0 


1.8 


6.2 


4.8(1.0) 


Cha3-C25 


25 


7.2 x 2.1 


3.4 


5.4 


4.3 


0.032 


0.081 


0.052 


62(17) 


0.18 


2.9 


4.9 


3.5 


8.1(2.2) 


Cha3-C26 


26 


2.3 x 2.1 


1.1 


5.2 


3.7 


0.031 


0.039 


0.054 


58(16) 


0.16 


2.8 


12.9 


3.6 


7.6(2.1) 


Cha3-C27 


27 


6.7 x 2.1 


3.2 


5.1 


2.3 


0.031 


0.071 


0.049 


63(18) 


0.17 


2.7 


4.9 


3.3 


8.2(2.4) 


Cha3-C28 


28 


4.9 x 2.1 


2.3 


5.1 


3.2 


0.031 


0.056 


0.038 


81(27) 


0.16 


2.7 


6.2 


2.6 


10.7(3.6) 


Cha3-C29 


29 


5.3 x 3.0 


1.7 


5.1 


3.7 


0.031 


0.082 


0.061 


50(13) 


0.18 


2.7 


4.6 


4.2 


6.5(1.7) 



Notes. (fl) Numbering of Gaussclumps sources like in Table[2] ib) Deconvolved physical source size (FWHM) and aspect ratio (R d ) of the fitted 
Gaussian. The minimum size that can be measured is 2100 AU (see Sect. 14.3.21 . The aspect ratio is the ratio of the deconvolved sizes along the 
major and minor axes. {c) Peak H2 column density. The statistical rms uncertainty is 1.0 x 10 21 cnT 2 . {d> Visual extinction derived from 2MASS. 
{e) Mass in the central beam (HPBW = 21.2") (M pcilk ), total mass derived from the Gaussian fit (M tot ), and mass computed from the flux measured 
in an aperture of 50" in diameter (M SQ »). The statistical rms uncertainties of M pea k and M so » are 0.006 and 0.010 M Q , respectively. W) Mass 
concentration m pC ak/™ 50 ". The statistical rms uncertainty is given in parentheses. <s) Ratio M toi IM^, with Mbe the critical Bonnor-Ebert mass (see 
Sect. l4.3.6t . (/,) Beam-averaged free-particle density within the central beam (/i pC ak) and mean free-particle densities computed for the total mass 
("mean) and the mass M iQ n in the aperture of diameter 50"(« 50 "). The statistical rms uncertainties of « pcak and n 5Q /i are 5.4 x 10 4 and 6.9 x 10 3 
cnT 3 , respectively. (,) Density contrast « peak /" 50 " • The statistical rms uncertainty is given in parentheses. 



are about a factor of 2 lower than in Cha I (see Fig. 6e of Paper I 
for comparison). 



4.3.2. Sizes 

The source sizes along the major and minor axes before and 
after deconvolution are listed in Cols. 6 and 7 of Table [2] 
and Col. 3 of Table [3] respectively. Their distributions are 
shown in Figs. [6^ and b, respectively, along with the distri- 
bution of mean size (geometrical mean of major and minor 
sizes, i.e. ^jFWHM mi ^ x FWHM m i n ). The average major, mi- 
nor, and mean sizes are 62 "> 35 ", and 46 *H ", respec- 
tively. These angular sizes correspond to physical sizes of 9300, 
5300, and 6900 AU, respectively. Only two sources have a ma- 
jor FWHM size larger than 110", and no source has a mi- 
nor or mean FWHM size larger than 90". The results of the 
Monte-Carlo simulations of PapeiT in the elliptical case (see also 
Sect. 12.31 and Appendix lA.U imply that these sources, although 
nearly all of them are faint with a peak flux density lower than 
150 mJy/beam, are not significantly affected by the spatial filter- 



ing due to the sky noise removal, with less than 15% loss of peak 
flux density and size. 

Like for Cha I (see Paper I), the accuracy to which we 
can measure the size of a weak ~ 5cr unresolved source is 
4.2" (21.2/5). Therefore, faint sources with a size smaller than 
~ 25.4" cannot be reliably deconvolved and we artificially set 
their size to 25.4" to perform the deconvolution. As a result 
the minimum deconvolved FWHM size that we can measure is 
~ 14" (2100 AU). The average deconvolved mean FWHM size 
is 40 t\t ", i-e. 6000 AU (see Fig.|6p). It is only 8% smaller 
than the average deconvolved mean size of the Cha I sources so 
the conclusions of Paper I hold as well: the Cha III sources have 
similar physical sizes as the Perseus cores, are probably larger 
than the Serpens cores (maybe by a factor of 1.5-2), are cer- 
tainly larger than the Ophiuchus cores (by a factor of 2-3), but 
are significantly smaller than the Taurus cores (by a factor of 3). 

The comparison to the population of dense cores in the 
Pipe nebula is less straightforward because these cores were ex- 
tracted from an extinction map with Clumpfind and only the radii 
of their lowest contours are available, not their FWHM sizes 
(lAlves et al]l2"007t lLada et"aDl2008t iRathborne et alJl2009l) . In 
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Fig. 5. Extinction map of Fig. Q] (red contours) overlaid on the 
870 yum continuum emission map of Cha III (black contours). 
The contour levels of the extinction map start at 3 mag and in- 
crease by steps of 1.5 mag. The thicker red contours correspond 
to Ay = 4.5 and 7.5 mag. The contour levels of the 870 fim map 
are the same as in Fig. [2] plus a dotted blue contour at -a. The 
dotted line delimits the field mapped at 870 /im. The field of 
view of LABOCA and the angular resolution of the extinction 
map are shown in the upper right corner. 



addition, the resolution of the extinction map is 1', i.e. 7800 AU 
at a distance of 130 pc dLombardi et al.ll2006l). Still, s ince their 
mean radius is about 19000 AU ( Rathb orne et al.1 2009), it seems 
very likely that the Pipe dense cores traced by the extinction are 
larger by a factor of a few compared to the Cha III cores detected 
with LABOCA. 



4.3.3. Aspect ratios and orientations 

The distribution of aspect ratios computed with the deconvolved 
FWHM sizes is shown in Fig. [6};. Based on the Monte Carlo 
simulations of Paper I, we estimate that a faint source can re- 
liably be considered as intrinsically elongated when its aspect 
ratio is higher than 1 .4. 76% of the sources are above this thresh- 
old and can be considered as elongated, which is similar to Cha I. 
The average aspect ratio is 2.3 It is similar to the ones mea- 
sured in Cha I, Serpens, and Taurus, somewhat larger than in 
Perseus, and significantly larger than in Ophiuchus (see Table 7 
of Paper I). 

The distribution of position angles (Col. 8 of Table |2]i does 
not show any preferred direction. An inspection by eye does not 
reveal any particular alignment of the elongated sources with 



the putative filaments mentioned in Sect. l3.1l either, especially in 
field Cha3-Center (Fig. [3};). 



4.3.4. Column densities 

The median peak H2 column density of the starless sources in 
Cha III is 6.4x 10 21 crrT 2 (Fig.[6ji), nearly identical to the median 
value found in Cha I (6.5 x 10 21 cm"" 2 , see Paper I). However, 
the average peak H2 column density of the starless sources in 
Cha III (6.9 x 10 21 cm" 2 ) is 1.4 times lower than the average 
peak column density of the starless cores in Cha I. It is 5 times 
lower than in Perseus and Serpens, and 4 to 9 times lower than in 
Ophiuchus (see Table 7 of Paper I). It appears to be significantly 
lower than in Taurus too (by a factor of 3), but since the Taurus 
sample is not complete and the source extraction methods differ, 
we cannot draw any firm conclusion. 



4.3.5. Masses and densities 

The distribution of masses and free-particle densities listed in 
Table [3] are displayed in Fig. [7] The 5cr sensitivity limit used to 
extract sources with Gaussclumps corresponds to a peak mass 
of 0.030 M and a peak density of 2.7 x 10 5 cm" 3 , computed 
for a diameter of 21.2" (3200 AU). The median of the peak 
mass distribution is 0.039 M Q , implying a median peak den- 
sity of 3.4 x 10 5 cm" 3 (Figs. [7^ and e). We give the mass in- 
tegrated within an aperture of diameter 50" (7500 AU) in Col. 9 
of Table [3] which is the aperture used for Cha I (Paper I). This 
aperture is well adapted to the Cha III sample too, for the same 
three reasons: it nearly corresponds to the average mean, unde- 
convolved FWHM size (see Sect. 14.3.21 and Fig. |6^), it is not 
affected by the spatial filtering due to the sky noise removal (see 
Appendix A. 1 and Table A. 1 of Paper I), and it is still preserved 
in the sum map at scale 5 (see Appendix C and Table C.l of 
Paper I). The median of the mass integrated within this aper- 
ture is 0.09 M , corresponding to a median mean density of 
6 x 10 4 cm" 3 (Figs. |7p and f). The median values of the peak 
and aperture masses are nearly the same as those of the Cha I 
sample, but the mean values are smaller by about 30%. Within 
the uncertainties, the mass distributions of Cha I and III are very 
similar, the only main difference being that the four most mas- 
sive cores of Cha I have no counterpart in Cha III. 

Figure [7J; shows the distribution of total masses computed 
from the Gaussian fits (Col. 8 of Table [3}. The completeness 
limit at 90% is estimated from a peak flux detection threshold at 
6.3cr for the average size of the source sample {FWHM = 46"fl 
It corresponds to a total mass of 0.18 M , slightly lower than in 
Cha I (0.22 M , see Paper I). Like in Cha I, the median total mass 
is very similar (0.20 M ), which implies that only 50% of the de- 
tected sources are above the estimated 90% completeness limit. 
Th e mass completeness l imit is slightly better than that obtained 
bv iKonvves et al.l (12010) for their 11 deg 2 sensitive continuum 
survey of the Aquila Rift cloud complex (distance 260 pc) with 
Herschel, and abo ut a factor of 6 better than that obtained by 
lEnoch et"a"D (|2008) in Perseus (see Paper I for details). It is how- 
ever a factor ~ 20 worse than that obtained with Herschel for 
the Polaris flare re gion which is at the same distance as Cha III 
(lAndreetal.ll2010h . 



5 For a Gaussian distribution of mean value m and standard deviation 
cr, the relative population below m - 1.28cr represents 10%. Therefore 
our peak flux detection threshold at 5<r implies a 90% completeness 
limit at 6.3cr, with a the rms noise level in the sum map at scale 5. 
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The mean density of each source is estimated from its total 
mass as derived from the Gaussian fits and a radius set equal to 
^jFWHM m . di x FWHM mhl . It is given in Col. 13 of Tableland 
the distribution for the full Cha III sample is shown in Fig. |7g. 
The average and median mean densities are 4.2*M x 10 4 and 
3.3 x 10 4 cm -3 , respectively. For Chamaeleon I, the correspond- 
ing numbers are 4.8^6 x 10 4 and 3.6 x 10 4 cm~ 3 , respectively 
(not given in Paper I). The Cha I and III sources have thus very 
similar mean densities, a factor of ~ 5 higher than the mean 
densit y of the cores of t he Pipe nebula extracted from extinction 
maps (iLada et al.ll2008l) . 

We estimate the mass concentration of the Cha III sources 
from the ratio of the peak mass to the mass within an aperture 
of 50" (Col. 10 of Table [3]l which is relatively insensitive to the 
spatial filtering due to the data reduction. A similar property is 
the density contrast measured as the ratio of the peak density 
to the mean density within this aperture (Col. 15 of Table [3]). 
The statistical rms uncertainties on the peak mass and the mass 
within 50" are 0.006 and 0.010 M©, respectively, which means 
a relative uncertainty of up to 25% for the weakest source. The 
distributions of both ratios are shown in Figs. [7ji and h and their 
rms uncertainties are given in parentheses in Cols. 10 and 15 
of Table[3] The two outliers with the largest ratios are also those 
with the highest relative uncertainty (about 30%). The upper axis 
of Fig.[7jl, which can also be used for Fig.[7ji, displays the expo- 
nent of the density profile under the assumptions that the sources 
are spherically symmetric with a power-law density profile, i.e. 
p cc r~ p , and that the dust temperature is uniform. The median 
mass concentration and density contrast are 0.44 and 5.8, respec- 
tively, similar to the Cha I sample. This corresponds to p ~ 2.0, 
suggesting that most sources are significantly centrally-peaked 
(see Paper I for the caveats of this estimate). It is similar to the 
exponent of the singular isothermal sphere. 

The upper axis of Fig.|7Ji, which can also be used for Fig.[7Jl, 
deals with an alternate case where the sources have a constant 
density within a diameter Dfl at and a density decreasing as r~ 2 
outside, still with the assumption of a uniform temperature. 
Under these assumptions, the measurements are consistent with 
a flat inner region of diameter 16" at most (2400 AU) for a few 
sources, but most sources have Dfl at < 10" (1500 AU), or cannot 
be described with such a density profile. 

4.3.6. Mass versus size 

The distribution of total masses versus source sizes derived 
from the Gaussian fits is shown in Fig. [8^. About 50% of the 
sources are located between the 5cr detection limit (solid line) 
and the estimated 90% completeness limit (dashed line), sug- 
gesting that we most likely miss a significant number of sources 
with a low peak column density. Figure |8j) shows a similar di- 
agram for the deconvolved source size. If we assume that the 
deconvolved FWHM size is a good estimate of the external 
radius of each source, then we can compare this distribution 
to the critical Bonnor-Ebert mass that characterizes the limit 
above which the hydrostatic equilibrium of an isothermal sphere 
with thermal support only i s gravitational ly unstable. This rela- 
tion M BE (R) = 2.4 Ra 2 JG dBonnorll 19561) . with M BE {R) the to- 
tal mass, R the external radius, a s the sound speed, and G the 
gravitational constant, is drawn for a temperature of 12 K as 
a solid line in Fig. [8J5. We define a B E = M tot /M B E- Only one 

6 The relative uncertainty of the ratio is equal to the square root of 
the quadratic sum of the relative uncertainties of its two terms, i.e. we 
assume both terms are uncorrelated. 



source (Cha3-Cl) has <?be > 1> i- e - is located above the criti- 
cal mass limit (see Fig. [8J5 and Col. 11 of Table |3). The thresh- 
old c^be = 0.5 approximately defines the limit above which an 
isothermal sphere is gravitationally bound if it is only supported 
by thermal pressure and the confinement by the external pres- 
sure is negligible. Only one source in addition to Cha3-Cl falls 
above this limit and may be gravitationally bound (Cha3-C2, 
see dash-dotted line in Fig. [8J3). Most sources, however, have 
a mass lower than the critical Bonnor-Ebert mass by a factor of 
2 to 6. The uncertainty on the temperature (see Appendix B.2 of 
Paper I) does not influence these results much since, even in the 
unlikely case of the bulk of the mass being at a temperature of 
7 K, the measured masses would move upwards relative to the 
critical Bonnor-Ebert mass limit by a factor of 1 .9 only, because 
the latter is also temperature dependent. 

The critical Bonnor-Ebert mass can also be estimated from 
the ex ternal pressure with M BE (P ext) = 1.18a 4 G~2 p^ t dBonnorl 
1 19561) . the external pressure being estimated from the extinc- 
tion of the environment in which the sources are embedded (see 
Paper I for the equations and references). The single source with 
a mass larger than M B E(P ext ) is the same as for M B e{R). The 
agreement between both estimates of M B e suggests that our esti- 
mates of the external radius and external pressure are consistent. 
No additional source falls above the threshold a B E - 0.5 based 
on M B E(Pext)- In summary, only one source is likely above the 
critical Bonnor-Ebert mass limit (Cha3-Cl), and one additional 
source may be gravitationally bound if it is supported by thermal 
pressure only (Cha3-C2). The implications of this analysis will 
be discussed in Sect. [5] 

The mass concentration Cm is plotted versus source size in 
Fig. [8}:. Cm is actually equal to the ratio of the peak flux to the 
flux integrated within the aperture of diameter 50". When the 
sources do not overlap, this ratio is nearly independent of the 
Gaussian fitting since the second and third stiffness parameters 
of Gaussclumps were set to 1, i.e. Gaussclumps was biased to 
keep the fitted peak amplitude close to the observed one and 
the fitted center position close to the position of the observed 
peak. The dashed line shows the expected ratio if the (not decon- 
volved) sources were exactly Gaussian and circular and allows 
us to estimate the departure of the sources from being Gaussian 
within 50". Most sources have a mass concentration consistent 
with the Gaussian expectation, but many of them have a signif- 
icant uncertainty on Cm that prevents a more accurate analysis. 
The obvious outlier toward the lower left is source Cha3-C18, 
which has a strong neighbor significantly contaminating its flux 
within 50" (source Cha3-Cl). 

There is no obvious correlation between the total mass or 
FWHM size of the sources and the visual extinction of the en- 
vironment in which they are embedded (see Fig. |9j. A simi- 
lar conclusion was drawn for Cha I (Paper I) and for the five 
nearby molecular clouds Ophiuch us, Taurus, Perseus , Serpens, 
and Orion based on SCUBA data (ISadavov et al.ll2010h . 

4.3.7. Core mass distribution (CMD) 

The mass distribution of the 29 starless sources is shown in 
Fig. [10] Its shape looks very similar to the shape of the mass 
distribution found in other star forming regions with a power- 
law-like behavior at the high-mass end and a flattening toward 
the low-mass end. In our case, the flattening occurs below the 
estimated 90% completeness limit (0.18 M ) and may not be 
significant. Above this limit, the distribution is consistent with a 
power-law but is very noisy. The exponent of the best power-law 
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Fig. 9. a Total mass versus visual extinction Ay for the 29 star- 
less sources found with Gaussclumps in the sum map of Cha III 
at scale 5. The dashed line shows the estimated 90% complete- 
ness limit (0.18 M ). b FWHM size versus visual extinction. 
The angular resolution (21.2") is marked by the dotted line. 
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Fig. 10. Mass distribution dA/dlog(M) of the 29 starless 
sources. The error bars represent the Poisson noise (in VA). The 
vertical dotted line is the estimated 90% completeness limit. The 
thick solid line is the best power-law fit performed on the mass 
bins above the completeness limit. The best fit exponent a\ og is 
given in the up per right corn er. The IMF of single stars corrected 
for bi naries (Kroupa feOOlL K01) and the IMF of multiple sys- 
tems (IChabrieri ho05. C05) are shown in dashed (red) and dot- 
dashed (blue) lines, respectively. They are both vertically shifted 
to the same number at 2 M . The do tted (purple) curve is the 
typica l mass spectrum of CO clumps (lBhtzlll993l : lKramer et alJ 
1998). 



fit (a = -3.0+0.8 for dA/dM,a log = -2.0±0.8 for dA/ d log(M)) 
is con sistent, within the uncertainty, with the value of ISalpeten 
(fl95l that characterizes the high-mass end of the stellar initial 
mass function (a = -2.35). However, it is also consistent within 
2cr with the exponent of the typical mass spectru m of CO clumps 
(a = -1.6, see lBlit3[l993l: iKramer et alj|1998l) . The sample is 
too small to distinguish statistically between these two types of 
mass distribution. 



4.4. Spatial distribution 

The distribution of nearest-neighbor projected distance is pre- 
sented in Fig.QT]for both samples of starless sources in Cha III 
and I. The median distance d m is 0.11 pc in Ch a III, a factor 
of 2 l arger than in Cha I (0.056 pc). We follow Gome zeit alJ 
(1993) to estimate the corresponding distribution for a sample 
of sources that would be randomly distributed in the plane of 
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Fig. 11. Distribution of nearest-neighbor projected distance for 
the starless sources of Cha III (a) and Cha I (b). In each panel, 
the dashed histogram shows the distribution expected for the 
same number of objects randomly distributed in the same area. 
The median value of each distribution is marked with a dotted 
line. 



the sky over the same area, assumed to be the surface of a disk 
of diameter equal to the largest projected distance between two 
sources (4.9 pc for Cha III and 4.0 pc for Cha I). The median 
distances of these random distributions are 0.38 pc for Cha III 
and 0.21 pc for Cha I, i.e. nearly a factor of 4 larger than ob- 
served in both cases. The starless sources in both clouds are thus 
significantly clustered. Assuming that the nearest-neighbor pairs 
are randomly oriented in the three-di mensional space, the true 
separation in three dimensions is -d m dGomez et alj|1993h . i.e. 
0.14 pc for Cha III and 0.07 pc for Cha I. 

Given that the median extinction of the local ambient 
medium in which the starless sources in Cha III are embedded is 
a factor of 2 lower than in Cha I (see Sect. 14.3. U . the density of 
this local ambient medium is most likely also lower, maybe up to 
a factor of 2. The characteristic length of thermal fragmentation 
is inversely proportional to the square root of the density. The 
difference in median nearest-neighbor separation seen between 
Cha III and I is therefore somewhat more pronounced than what 
could naively be expected from thermal fragmentation. 

The remarkable alignment of six 870 fim sources in field 
Cha3-Center (see Fig. [3};) deserves some further analysis. The 
sources are nearly uniformly distributed along a straight line, 
with a mean projected separation of 0.16 ± 0.03 pc. For a 
density of the intercore medium of 4 x 10 3 cirT 3 in this 
region (iGahm et alJ 120021) . the Jeans length 2na s / yjAnGp is 
about 0.36 pc at 12 K. The intercore separation would match 
this Jeans length if the inclination angle of the putative fila- 
ment was 26°, which is statistically unlikely. However, the ef- 
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fective scale of fragmentation in a magnetized and/or rotat- 
ing filament is expected to decrease wi th increasing magnetic 
field strength and/or ro tation level (e.g. INakamura et al. 1993; 
Matsumot o et alll994l) . The measured core separation could be 
used to estimate the magnetic field strengt h and/or rotation leve l 
as was done for a filament in Orion A (lHanawa et al.l IT993b . 
With 

^fwhm the full width at half maximum of the filament, 
A the projected core separation, a and £1 the ratios of the mag- 
netic pressure and centrifugal force to the thermal pressure, re- 
spectively, and i the incl ination of the filame nt along the line 
of sight, Equation 13 of lHanawa et al.l d 1993b yields a + ft > 



(1.75 x 



+ 0.6) - 10. The equality holds for a pitch 



4.3 ,1/sini 

angle 8 = 0°, with 9 characterizing the relative strength of the 
poloidal to axial magnetic fields. The width of the filament can- 
not be estimated from t he LABOCA map . With the C ls O 1-0 
map shown in Fig. 9 of Gah m et al.l (120021) . we roughly estimate 
4whm ~ 0.25-0.40 pc, which is larger than the typical width 
of the filaments recently detected with Her schel in three other 
nearb y clouds (median width 0. 1 pc, see Arzoumanian et al. 
1201 lb . Since our estimate from C 18 is rather uncertain, we 
also consider below (in parentheses) the case where the width 
of the filament is 0.1 pc. If the filament is in the plane of the 
sky (z = 90°) and 6 = 0°, we obtain a + /3 ~ 53-180 (6 with 
0.1 pc). Since t here is apparently no significant level of rotation 
in this filament dGahm et al.ll2002l) , we derive a ~ 53-180 (6). If 
the filament is not in the plane of the sky, a would be less than 
1 for i $ 7.5-12° (31°). An inclination to the line of sight of 
7.5-12° is statistically unlikely, so we are left with a very high 
level of magnetic pressure or an overestimated filament width. 
If this alignment of 6 regularly spaced cores is really the result 
of fragmentation in a filament, then we conclude that either the 
filament is strongly magnetized, or it is much thinne r than it ap- 
pears in C ls O 1-0, or the model of INakamura et al .1 (119931) of a 
magnetized, self-gravitating, isothermal filament in equilibrium 
does not apply to that filament. Alternatively, the observed reg- 
ular structure may have nothing to do with fragmentation and 
simply represent transient periodic overdensities produced by 
gravitational-m agnetoacoustic waves that will be damped away 
(lLangerlll978l) . 



5. Discussion 

5.1. A puzzling population of starless cores in Cha I and III 

Based on the comparison to the Bonnor-Ebert mass limit, we 
estimate that only one (or at most two) source(s) out of 29 is a 
candidate prestellar core in Cha III (see Sect. l4.3.6i l. This yields a 
fraction of candidate prestellar cores of 3-7%, a factor of 2 lower 
than in Cha I (5-17%, see Paper I). Apart from the few candi- 
date prestellar cores, the population of starless cores in Cha III is 
very similar to the one in Cha I since they have nearly the same 
median peak, aperture, and total masses as well as the same me- 
dian size and aspect ratio (compare Figs. [6] and [7] to Figs. 7 and 
8 of Paper I). They also follow the same correlation in an q-be 
versus M tot diagram (see Fig. ITSb . The main striking difference 
is that the visual extinction of the medium in which the Cha III 
sources are embedded is on average a factor of 2 lower than in 
Cha I. Although we a priori cannot exclude that the extinction 
laws of both clouds may differ or that there may be more con- 
tamination by foreground stars toward Cha III leading to an un- 
derestimate of the extinction, we rather consider that this differ- 
ence may come from the density structure or the physical pro- 
cesses at work in the clouds. As mentioned in Sect. [1] Cha III 
looks much more filamentary in cold dust emission at 100 fim 
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Fig. 12. Ratio of total mass to Bonnor-Ebert critical mass as a 
function of total mass. The filled diamonds and open squares 
show the Cha I and III sources, respectively. The dashed line 
shows the limit above which a source is gravitationally unstable. 
The gravitationally bound sources are located above the dotted 
line, provided they are supported by thermal pressure only and 
the external pressure is negligible. 



than Cha I (see Fig. 7 of Boulan ger et al.| [l998). If this is also 
the case on scales smaller than the resolution of our extinction 
maps (3'), one could expect lower extinctions for the "ambient" 
medium in Cha III compared to Cha I. Alternatively, this differ- 
ence in ambient extinction (and thus ambient density) for two 
otherwise similar populations of starless cores may indicate that 
the structuring of the interstellar medium into seeds of cores does 
not depend much on the local gravity and may be dominated by 
other processes su ch as turbulence and magnetic fields (see also 
lAndr6 et al.|[20lob . 

The sample of starless sources in Cha III being small, the 
accuracy of the CMD is not sufficient to compare it with the 
IMF and the CO clump mass spectrum (see Sect. I4.3.7l i. To en- 
large the statistics, the source samples of Cha I and III can be 
merged, both surveys having nearly the same mass complete- 
ness limit and their populations of starless cores having very 
similar properties. The full sample contains 89 sources (60 in 
Cha I and 29 in Cha III) and its mass distribution is shown in 
Fig. [T3h . We can also consider only the starless sources of this 
enlarged sample that are not candidate prestellar cores, i.e. those 
with ffBE < 1 ■ This yields a sample of 85 sources (57 in Cha I 
and 28 in Cha III), the excluded sources being Chal-Cl-3 and 
Cha3-Cl . The mass distribution of this slightly smaller sample is 
presented in Fig. [TJb . Above the 90% completeness limit, both 
distributions are well fitted by a single power-law with an ex- 
ponent ai g = -1.5 ± 0.3 and -1.7 ± 0.4, respectively. This is 
steeper than, but still consistent within lcr with, the Salpeter ex- 
ponent of the high-mass end of the stellar IMF, and it is defi- 
nitely much steeper than the CO clump mass distribution that is 
characterized by a\ og — -0.6 CBhtz1 ll993l : iKramer et al. 1998fl 
Even if the combined sample of Cha I and III starless sources 



7 We note that llkeda & Kitamural feOllI) fi nd a = -2.1 ± 0.2 with 
C ls O 1-0 in S140, the same cloud for which IKramer et alj dl998h de- 
rived a = — 1 65 ± 0.18 with C 18 2-1. They argue that the study of 
Kram er et al] (1998) is limited to the central part of the cloud and is 
likely biased toward high-mass cores (hence yielding a flatter CMD). 
They also show that a poor spatial resolution leads to underestimating 
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a log = -1.5 ± 0.3 




a log = -1.7 ± 0.4 
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Fig. 13. Mass distribution dAVdlog(M) of the starless sources 
of Cha I and III. Panel a shows all sources while panel b dis- 
plays only those with M/M^e < 1 -0. The error bars represent 
the Poisson noise (in yN). The vertical dotted line is the av- 
erage (0.20 M Q ) of the estimated 90% completeness limits for 
Cha I and III. The thick solid line is the best power-law fit per- 
formed on the mass bins above the completeness limit. See cap- 
tion of Fig.[T0]for other details. The K01 and C05 IMFs are both 
vertically shifted to the same number at 5 M . 



detected with LABOCA is large (89 sources), there are only 
about 50 sources above the completeness limit. We thus expect 
that the ongoing sensitive survey of the Chamaeleon clouds per- 
formed with Hers chel in the frame of the Gould Belt Survey 
dAndre et al.l20 10) will provide a more complete sample of star- 
less cores and thus a more robust CMD. In particular, the shape 
of the mass distribution of the starless sources that are not can- 
didate prestellar cores will be of prime importance. 

The population of starless cores detected with LABOCA in 
Cha I and III is very puzzling: although most of these sources do 
not appear to be prestellar based on the Bonnor-Ebert criterion, 
their mass distribution seems to be consistent with the stellar 
IMF at the high-mass end. For Cha I, we argued that a loss of 
thermal support via cooling would not be sufficient to bring these 
cores above the Bonnor-Ebert mass limit, hence that they are 
unlikely to form stars in the future (Paper I). Along with other 
arguments, this suggested that star formation is over in Cha I. 
In the same way, we could a priori conclude (but see the next 
sections) that, apart from the northern part of the cloud where 
one candidate prestellar core is found, Cha III does not seem to 
be able to start the process of star formation. 



\a\ if the resolution is worse than 0.1 pc. The "classical" index a = -1.6 
is therefore most likely valid for larger CO clumps only. 
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The level of turbulence is very similar in both clouds, as mea- 
sured via the linewidths of the 7 = 1-0 transitions of CO and 
C ls O (see Sect. l5.3l below). Taken at face value, turbulence thus 
does not seem to be the key parameter promoting star formation 
since Cha I formed many stars in the past while Cha III did not. 
However, the present level of turbulence in Cha I may not reflect 
the initial conditions when the first stars were formed. On the 
one hand, the level of turbulence in this cloud may have been 
lower in the past and have been raised as a result of stellar feed- 
back via, e.g., molecular outflows, possibly preventing the pro- 
cess of star formation to continue at the present epoch. On the 
other hand, if the feedback of the Cha I (low-mass) YSOs has 
not been sufficient, the turbulence may have decayed with time. 
In this case, the same behaviour would a priori be expected for 
both clouds, which would not explain their difference in terms 
of past star formation activity. 

5.2. Are the Chamaeleon cores similar to the Pipe cores? 

At this stage, it is very instructive to compare the population 
of starless cores in Cha I a nd III to the one found in the Pipe 
nebula via extinction maps dAlves et al. 2007t lLadaet"aiT 2008: 
Rathborne et al. 2009). The Pipe nebula is the only other nearby 
star forming region in the pre-Herschel era for which a popula- 
tion of dense cores that are mostly gravitationally unbound has 
been found, with a mass distribution still consistent with the stel- 
lar IMF at the high-mass end. In addition, its CMD departs from 
a single power law and has a break like the IMF but at a higher 
mass of 2.7 + 1.3 M a, suggesting a star for mation efficiency of 
~ 20% in this cloud (Rathbo rne et al.l l2009). Interestingly, only 
the sources with a mass higher than this characteristic mass ap- 
pear to b e prestellar based on the Bon nor-Ebert criterion (see 
Fig. 9 of lLada et al.ll2008l) . lLada et al l (120081) suggest that the 
gravitationally unbound cores are pressure confined, the exter- 
nal pressure being most likely provided by the weight of the 
cloud itself . They thus do not appear to be transient structures. 
Lada et al. I J2008h consider two mechanisms that may turn these 
stable but unbound cores into prestellar cores: either an increase 
of the external pressure produced by the contraction of the whole 
cloud, or a mass increase if these cores have not obtained their fi- 
nal mass yet. The former mechanism would require a decrease of 
a factor of 2 in cloud radius, which may be possible via the dissi- 
pation of its supersonic turbulence but casts a potential timescale 
issue since the more massive unstable cores would form stars 
much more rapidly than the less massive ones (see [Lada et al] 
2008 for details). 

The starless cores in Cha I and III that are presently unstable 
according to the Bonnor-Ebert criterion are also the most mas- 
sive ones, like in the Pipe nebula, but the transition occurs at a 
lower mass (~ 1 versus ~ 3 M , see Fig. [T2l . With a median 
mass of 0. 9 M H and a median rad ius of 0.07 pc (computed from 
Table 2 of Rat hborne et al.lr2009l) . the Pipe cores are bigger than 
the Chamaeleon ones by a factor of ~ 4 in mass and ~ 2.5 in 
radius (but see the caveats about the radius in Sect. 14.3.21 . While 
this could be an intrinsic property, we believe that it results from 
the different tracers used to extract the cores: extinction maps are 
more sensitive to extended, low-density material than 870 /-im 
dust emission maps. The factor of 3 difference between the tran- 
sition masses mentioned above may therefore simply result from 
the fact that the LABOCA masses do not include the low-density 
material surrounding each core. 

A major difference exists however between the CMDs of the 
Pipe nebula and the Chamaeleon I and III clouds: there is no ev- 
idence for a break in the Chamaeleon CMD down to the 90% 
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completeness limit of 0.2 M , while such a break is seen around 
2.7 M Q in the Pipe nebula. Even the possible mass scaling factor 
of 4 mentioned above would not be sufficient to explain this dif- 
ference. However, the low-density material possibly missing in 
the LABOCA masses could affect the shape of the CMD rather 
than contributing as just a scaling factor. On the one hand, if the 
Chamaeleon cores are seeds of prestellar cores in the process of 
accumulating mass (see Sect. l5.3l below). we may conjecture that 
the break in the CMD could arise from this mass accumulation 
process. On the other hand, a large fraction of the lowest mass 
cores may never become prestellar and they may currently hide 
the true shape of the prestellar CMD in Chamaeleon. 



5.3. Can the starless cores in Cha I and III turn prestellar? 

The possibility that the unbound starless cores of the Pipe and 
Cha I/III clouds are still accumulating mass has been mentioned 
in Sect. l5.2l and is very attractive in light of the results obtained 
by recent numerical simulations. In this section, we compare the 
properties of Cha I and III to several kinds of numerical simula- 
tion. 

iGomez et al.l (12007) study the formation and collapse of qui- 
escent cloud cores induced by focused compressions (or "con- 
vergent flows") in a cloud of diameter 1 pc that initially has 
a constant sub-Jeans density of 113 cm -3 and a uniform tem- 
perature of 1 1 .4 K. The velocity amplitude of the compressing 
impulse is 0.4 km s _1 (Mach number of 2). A mild shock prop- 
agates inwards, the material left behind it (the envelope) being 
set into infall motions. The shock bounces off the center and ex- 
pands outwards, leaving a quiescent core in the inner part. The 
structure of the core+envelope system then resembles a trun- 
cated Bonnor-Ebert sphere with a flat inner density profile and a 
density falling as ~ r~ 2 in the (infalling) envelope. Depending on 
the initial conditions (position of the impulse within the cloud), 
the core can gain enough mass from the envelope to become 
gravitationally unstable and start collapsing. In this model, the 
inner quiescent core with a density of ~ 10 5 cm -3 is initially un- 
bound and confined by the ram pressure of the inflowing gas. 
It grows in size and mass until it gets dominated by gravity. 
Interestingly, there is a significant time delay of ~ 5 x 10 5 yr 
between the formation of the inner core (when the system starts 
to look like a pseudo Bonnor-Ebert sphere) and the onset of grav- 
itational collapse, due to the growth in mass. 

The C ls O 1-0 linewidths in Cha I and III ar e on the order 
of 0.8 -0.9 km s _1 at an angular resolution of 2.7' dMizuno et al.l 
Il999h . implying an rms (turbulent) velocity dispersion of 0.3- 
0.4 km s . The CO 1-0 linewidths, tracing ev en lower-density 
mater i al, are a factor of ~ 2.5 times larger (iBoulanger et al.l 
1998t Mizuno etal. 120011) . The initial conditions chosen by 



Gomez et al.l ([2007) are therefore plausible for Chamaeleon 



from a kinematic point of view. The typical densities of the 
Chamaeleon starless cores are also similar to the density of the 
inner core formed in these simulations. If this scenario holds for 
Chamaeleon, a fraction of these presently unbound cores could 
turn prestellar in the future (in less than ~ 5 x 10 5 yr) if they 
gain enough mass to become unstable. We note, however, that 
there is no evidence for significantly flattened inner density pro- 
files in the Chamaeleon sample (see Sect. 14. 3.31 and Paper I), but 
this may not be a shortcoming for this sce nario, the spherically 
symmetric model and initial conditions of Gom ez et al.l (2007) 
being highly i dealis ed. Finally, since some of the simulations of 
Gomez et al. 1 (120071) do not lead to the formation of an unsta- 
ble core prone to collapse, a fraction of the Chamaeleon starless 



cores may never reach the critical mass and simply get dispersed 
in the future (see their simulation SI). 

Another fru it ful a pproach is the one followed by 
Clark & Bonne H 42005). These authors examine the for- 
mation of bound coherent cores in molecular clouds supported 
by (decaying) turbulence. They model a cloud of mass 32.6 M Q 
and diameter ~ 0.3 pc at 10 K with an initially uniform density 
of 5.6 x 10 4 crrT 3 and turbulent velocities characterized by 
an initial effective Mach number of 5.3 and a power spectrum 
P(k) oc AT 4 (their simulation 2). A large number of fragments^ 
are formed, most of them being initially unbound. Their mass 
distribution is well fitted by a Salpeter slope at the high-mass 
end when star formation sets in (at t ~ tg, with tg the initial 
free-fall time of the cloud). Interestingly, the mass distribution 
of these mostly unbound fragments is slightly steeper at the 
high-mass end at an e arlier stage (t - 0.6%, see Fig. 3 of 
IClark & Bonnelllr2005l) . A second important point is that only 
the most massive fragments are gravitationally bound (see their 
Fig. 5). In this simulation, the unbound fragments grow in mass, 
as a result of coagulation, but only a small fraction of them 
become gravitationally unstable and can start collapsing. The 
number of gravitationally unstable cores formed in this way (11 
for their simulation 2) is on the same order of magnitude as the 
initial number of mean Jeans masses in the cloud (33). 

The population of starless cores in Cha I a nd III share sim- 
ilar pr operties with the simulated fragments of IClark & Bonnell 
(2005]): their CMD resembles the IMF at the high-mass end, 
most cores are unbound, and only the most massive ones are 
gravitationally unstable. The slope of the Chamaeleon CMD at 
the high-mass end is even slightly steeper (but only at the l<x 
level) than the Salpeter value, like for the simulated fragments 
before star formation sets in. However, most objects qualified as 
fragments in the simulation are much too small and too faint 
to be detected in our continuum maps. There are more than 
1000 fragments formed in the simulated cloud of projected area 
0.053 pc 2 , while we detect only 89 cores in 17 pc 2 . Because 
of the limited sensitivity and angular resolution of our maps, the 
low number of detected starless cores does not really rule out this 
scenario of turbulence-generated, small, unbound fragments that 
coagulate with time. We may be seeing only those fragments that 
have already grown enough by coagulation to be detected. One 
issue with this analogy may b e that the initial level of t urbulence 
assumed in the simulations o flclark & Bonnel (l2005h is signif- 
icantly higher than the one characterizing the regions of Cha I 
and III that have similar densiti es of ~ 5 x 10 4 cm' 3 (Mach num- 
ber < 2 based on C 18 1-0, see lMizuno et al.lll999tlGahm et al.l 
120021) . 

The comparison of the properties of the Cha maeleon starless 
cores to th ose produced in the simulations of IClark & Bonnelll 
(2005]) and IGomez et al.l d2007l) leaves open the hypothesis that 
a fraction of these starless cores may become unstable in the 
future by accumulating mass from larger scales. Since most 
Chamaeleon starless cores are located in filamentary structures, 
the next question is whether this gain in mass could simply oc- 
cur along the filaments or would require accumulation of ma- 
terial from larger scales. In Cha I, the mass of the filaments 
traced with LABOCA is, without counting the cores embedded 
within them, roughly equal to the mass of the starless cores. If 
this LABOCA mass is representative of the real mass of the fil- 



IClark & Bon nell (2005) use the word "clumps" for these objects, 
but we prefer to call them "fragments" since, in the context of observa- 
tions, the word "clump" is usually gi ven to larger-scale structures that 
may contain several cores (see, e.g.. lWilliams et al.ll2000h . 



13 



A. Belloche et al.: Will the starless cores in Chamaeleon I and III turn prestellar?' 



aments, then the mass gain could not be larger than a factor of 
2 and would not be sufficient to bring most starless cores above 
the Bonnor-Ebert limit, provided the product of their temper- 
ature and radius does not at the same time decrease by more 
than a factor of 2 (see the definition of M^(R) in Sect. 14.3.6b . 
In their simulations of d ense core formation in supersonic tur- 
bulent converging flows, iGong & Ostrikerl (1201 ll) report that fil- 
aments are formed in post-shock regions at the same time as 
overdense regions within these filaments condense into cores. It 
is thus tempting to conclude that the starless cores in Cha I and 
III may grow in mass at the same time as the filaments do un- 
der the influen ce of turbulence and self- gravity. Interestingly, the 
simulations o f lGong & Ostrikerl (1201 lh can produce alignments 
of nearly regularly spaced cores. For a Mach number of 5, the 
core separation is on the order of 0.2 code units, i.e. 0. 15-0.47 pc 
for a mean density of free particles of 4000^-00 cirT 3 (see their 
Fig. 7, top left and bottom right panels). This is very much remi- 
niscent of the remarkable alignment of regularly spaced cores in 
Cha3-Center (see Sect. l4.4l andFig.[3b). The CO linewidths mea- 
sured on large scales in Ch a III are on the order of 2-2.5 km s' 1 , 
i.e. a Mach number of 4-5 dBoulanger et al.l 1998; Mi zuno et alj 
1200 ll) . Depending on the scale and the assumptions made, the es- 
timates of the mean density range bet ween 400 and 400 cm -3 
for Cha3-Center (see Sect. EHl and iGahm et al.l l2002h . Given 
the unknown inclination of the filamentary structure, the pro- 
jected separation of the aligned cores in Cha3-Center could be 
in rough agreement with the one produced in the simulations of 
Gong & Ostrikei] d201 ll) . Since some of the cores become gravi- 
tationally unstable in these simulations, this may also be the fate 
of some of the starless cores in Cha3-Center, and by extension 
in the full sample of Chamaeleon starless sources. 

5.4. A large fraction of "failed" cores in Cha I and III? 

In the simulations of [Gomez et aTJ (120071) . the phase of mass 
growth, from the time of central core formation on, lasts about 
as long as the collapse phase. In their ID sp herically symmet- 
ric sim ulations of converging supersonic flows. lGong & Ostrikerl 
(2009) find a duration of the core-building phase about 9 times 
longer than the duration of the colla pse phase, which is als o 
confirmed by their 3D simulations (Gong & Ostrikerl 1201 ll) . 
How ever only a fraction of t his core-building phase is observ- 
able. |Gong^&i)!trikel (|200^) estimate that 10-20% of the core 
building phase is observable if one takes as observability cri- 
terion a density contrast between the center and the edge of a 
core, Pc/Pedge, higher than 5. This means a duration of the ob- 
servable core-building phase 1-2 times longer than the duration 
of the collapse phase, which is rougly consistent with the result 
of Gomez et al. f (l2007l) . 

If we associate the foui0 starless cores of Cha I and III that 
have c^be > 1 with the collapse phase and the other starless cores 
(84) with the observable core-building phase, then we would ex- 
pect only 4-8 additional cores (3-6 in Cha I and 1-2 in Cha III) 
to become gravitationally unstable in the future, provided the 
process of star formation occurs at a constant rate and the model 
prediction of the duration of the core-building phase mentioned 
in the previous paragraph is correct. If the LABOCA masses are 
underestimated by a factor of 2, then there are about 10 and 2 
candidate collapsing cores at the present time, plus 10-20 and 
2-4 additional starless cores that could become gravitationally 
unstable in Cha I and III, respectively. Given these numbers, at 



9 We exclude Chal-Cl, i.e. the candidate first hydrostatic core Cha- 
MMS1, but we include Cha 1-C4 that has qtbe > 1 based on the pressure. 



least 29-39 and 23-25 starless cores, i.e. 49-66% and 79-86% 
of the population of starless cores detected with LABOCA, will 
most likely never become gravitationally unstable in Cha I and 
III, respectively. These fractions may even be underestimated if 
the cores detected with LABOCA have a ratio p c /p e dge signif- 
icantly higher than 5, which may well be the case since c„ in 
Table[3](and Table 6 of Paper I) is likely smaller than the true ra- 
tio Pc/pedge- However, these large fractions of "failed" cores hold 
only if the simulations of converging flows mentioned above are 
appropriate descriptions of the physical processes at work in the 
Chamaeleon molecular clouds. If the observable core-building 
phase lasts significantly longer than the collapse phase, then the 
actual fraction of "failed" cores could be smaller. 

In light of this new analysis, it cannot be excluded that Cha I 
will go on forming stars. However, even if mass growth occurs 
and the number of prestellar cores increases by a factor 2-3 as 
suggested in the previous paragraphs, this number will still be 
too low by a factor 2-4 (if the LABOCA masses are correctly 
estimated) to be consistent with a constant star formation rate 
since the time when the current pre-main-sequence stars were 
formed (2 Myr ago, see Paper I and references therein). In ad- 
dition, the lack of Class protostars by a factor maybe as large 
as 10 in Cha I (see Paper I) remains a strong indication that the 
star formation rate has decreased with time in this cloud. Cha I 
may not be at the end of the process of star formation as was 
suggested in Paper I, but it is at least secure to conclude that its 
star formation rate has decreased with time by a factor > 2-4 
over the last 2 Myr. On top of this overall decrease, it may have 
significantly fluctuated, by a factor 2-5 to account for the lack 
of Class protostars. 

5.5. A prime target to study the formation of prestellar cores 

With one candidate prestellar core and one or two additional star- 
less cores that could turn prestellar (see Sect. 15.41 ). Cha III seems 
to be able to form stars, even if presently at a very low rate. This 
is somewhat in contrast with the Polaris flare region, another 
presently non-star-forming molecular cloud. A large population 
of starless cores was recently uncovered with Herschel in that 
cloud but the cores typically lie one to two orders o f magnitude 
below the Bonnor-Ebert mass limit (see Fig. 4 of lAndre et al.l 
120101) . Since the bulk of the population of starless cores detected 
with LABOCA in Cha III lies a factor 3-4 below the Bonnor- 
Ebert mass limit, the tip of the Polaris population has otbe lower 
by a factor ~ 3 compared to the Cha III starless cores. This sug- 
gests that Cha III may be closer to fo rm prestellar cores than 
Polaris (but see Heitha usen et al.l l2002 for possible evidence of 
gravitationally bound cores in the latter). 

With a large fraction of starless cores presently not prestellar 
but still some evidence that star formation can occur, Cha III be- 
comes a prime target to study the formation of prestellar cores, 
i.e. the core-building phase when mass is accumulated, and thus 
the onset of star formation. In particular, it will be essential to 
get observational constraints on the duration of the core-building 
phase prior to the phase of gravitational collapse. In the previous 
section, we relied on predictions of numerical simulations to es- 
timate this duration and derive the fraction of "failed" cores, but 
this needs to be tested observationally. 

Finally, on a very speculative level, the factor of 2 differ- 
ence between the median nearest-neighbor distances of the star- 
less cores in Cha I and III (0.07 pc versus 0.14 pc, respectively, 
see Sect. 14.41 and Fig. [TTJ could be interpreted as fragmenta- 
tion being more efficient or more advanced in Cha I compared 
to Cha III, which could in turn have played a role in the fact that 
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Cha I has been more efficient than Cha III in forming stars. In 
that respect, it would be instructive to investigate whether signif- 
icant differences in, e.g., magnetic field structure exist between 
both clouds that could explain this speculative difference in their 
ability to fragment. 

6. Conclusions 

We performed a deep, unbiased, 870 /im dust continuum sur- 
vey for starless and protostellar cores in Chamaeleon III with the 
bolometer array LABOCA at APEX. The resulting 0.9 deg 2 map 
was compared with a map of dust extinction. The analysis was 
performed by carefully taking into account the spatial filtering 
properties of the data reduction process following the prescrip- 
tions of Paper I. The extracted sources were compared to those 
found in Cha I (Paper I) and other molecular clouds. Our main 
results and conclusions are the following: 

1 . The mass detected with LABOCA (23 M ) represents only 
6% of the cloud mass traced by the dust extinction, and about 
54% of the mass traced by the C ls O 1-0 emission. 

2. 29 sources were extracted from the 870 /im map, all of them 
being starless. No unresolved source is detected, which is 
consistent with the absence of any known young stellar ob- 
ject in this cloud. 

3. The starless sources are found down to a visual extinction of 
1.9 mag. Unexpectedly, about half of the sources are found 
below 5 mag, which is in marked contrast with other molec- 
ular clouds, including Cha I, where starless cores are only 
found above this threshold. Since the LABOCA surveys to- 
ward Cha I and III have the same sensitivity and were anal- 
ysed in the same way, this result may point to an intrinsic 
structural difference between the two clouds. 

4. The 90% completeness limit of our 870 /im starless core sur- 
vey is 0.18 Mq. Only 50% of the detected starless cores are 
above this limit, suggesting that we may miss a significant 
fraction of the existing starless cores. 

5. Although the distribution of starless sources suggests the ex- 
istence of filaments, these filaments are not detected with 
LABOCA in Cha III, while the LABOCA map of Cha I 
shows clear evidence of filamentary structures with the same 
sensitivity. 

6. There is a remarkable alignment of 6 nearly equally spaced 
sources in Cha III which may have been produced by turbu- 
lent fragmentation of a filament or simply represent transient 
periodic overdensities. 

7. Apart from their distribution of ambient extinction, the 
Cha III starless cores share very similar properties with those 
found in Cha I. They are less dense than those in Perseus, 
Serpens, Ophiuchus, and Taurus by a factor of a few on av- 
erage. 

8. At most two sources (< 7%) are above the critical Bonnor- 
Ebert mass limit in Cha III, which suggests that a large frac- 
tion of the starless cores may not be prestellar. Only the 
most massive cores in Cha I and III turn out to be candidate 
prestellar cores according to the Bonnor-Ebert mass crite- 
rion, in agreement with the correlation observed in the Pipe 
nebula. 

9. The mass distribution of the 85 starless cores of Cha I and 
III that are not candidate prestellar cores is consistent with a 
single power law down to the 90% completeness limit, with 
an exponent close to the Salpeter value. There is no evidence 
for a break such as the one seen in the core mass distribution 
of the Pipe nebula, and in the stellar IMF. 



10. A fraction of the Cha I and III starless cores that are presently 
not candidate prestellar cores may still be growing in mass 
and turn prestellar in the future. Based on predictions of nu- 
merical simulations of turbulent molecular clouds concern- 
ing the duration of the core-building phase, we estimate that 
at most 50% and 20% of the starless cores detected with 
LABOCA in Cha I and III, respectively, may form stars. 

Given their large fraction of starless cores that do not appear 
to be prestellar yet, Cha I and III turn out to be excellent sites 
to study the formation of prestellar cores. Since we find some 
evidence that star formation can start in Cha III, this cloud even 
becomes a prime target to investigate the onset of star formation. 
The main next step to make further progress will be to get obser- 
vational constraints on the duration of the core-building phase 
that precedes the gravitational collapse. 
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Appendix A: LABOCA data reduction: spatial 
filtering and convergence 

A. 1. Spatial filtering due to the correlated noise removal 

In a previous work, we characterized in detail the spatial filter- 
ing of LABOCA data due to the correlated noise removal and 
the high-pass filter (Paper I). Since these observations were per- 
formed using a combination of rectangular and spiral scanning 
patterns, the latter being more compact (diameter 6'), the prop- 
erties of the spatial filtering affecting the Cha III dataset, that 
was obtained with rectangular scanning patterns only, could a 
priori differ: the correlated noise removal depends on the size 
of the camera only, but the baseline removal is performed scan- 
(or subscan-)wise, i.e. on scales that can be larger than the cam- 
era, so it depends on the scanning pattern. We performed Monte 
Carlo simulations to evaluate the filtering properties follow- 
ing exactly the same procedure as described in Appendix A of 
Paper I, with four sets of 25 circular, artificial, Gaussian sources 
each. Figures lA.lfa -e and Table I A. 1 1 present the results. As a 
control check, we show in Figs. IA.lt — i the results of the fits per- 
formed on the artificial sources directly inserted in the final un- 
smoothed continuum emission map of Cha III. 

To first order, the results listed in Table |A~T1 are very simi- 
lar to those obtained for Cha I. The difference between strong 
and weak sources concerning the peak flux density is less pro- 
nounced for Cha III than for Cha I but this may not be statis- 
tically significant because the estimate for the weak artificial 
sources in Cha III is based on very few extended sources: there 
are only 3 such sources between 130" and 260", two of them 
behaving like the sample of stronger sources and one being sig- 
nificantly more affected by the filtering (see the outlier point at 
(145", 0.63) in Fig. IA.lb ). However, the size and integrated flux 
suffer more for Cha III than for Cha I. This is surprising since 
one would a priori expect the rectangular scanning pattern to 
be less sensitive to spatial filtering than the more compact spi- 
ral scanning pattern. The reasons for this difference are still un- 
clear. It may be due to a reduction of the number of well-working 
pixels in 2010 compared to 2007-2008: on average, the Cha III 
scans had 224 + 2 well-working pixels while the Cha I scans had 



Table A.l. Filtering of artificial Gaussian sources due to the cor- 
related noise removal depending on the input size (FWHM). 



Source sample 




Size where 




(1) 


ft, 

§■ = 85% 

Jp 

(2) 


-^-U = 85% 

S 

(3) 


f° 

f- = 85% 

'"" (4) 


Circular sources: 








strong sources 

(> 150 mJy/beam) 


230" 


210" 


140" 


weak sources 

(< 150 mJy/beam) 


230" 


110" 


60" 


all sources 


230" 


180" 


100" 



Notes. f° and f' p are the output and input peak flux densities, s' the 
input size (FWHM), s° and s° 2 the output sizes (FWHM), and f° m and 
f' mt the total output and input fluxes over an aperture of radius equal to 
the input FWHM. 




u 0.05 




10 15 
Iteration 



10 

Iteration 



Fig.A.2. a Convergence of the total 870 /im flux recovered in 
Cha III by the iterative process of data reduction, b Relative flux 
variation between iterations i and i — 1 as a function of iteration 
number i. 

232 + 18. Alternatively, this could suggest that mixing rectan- 
gular and spiral patterns is more robust to spatial filtering than 
scanning with rectangular patterns only. 

A.2. Convergence of the iterative data reduction 

Like for Cha I, the iterative process of the data reduction im- 
proves the recovery of extended emission. Figure lA~2l shows that 
the convergence of the total astronomical flux is reached within 
1% at iteration 20. Figure |A~31 shows how the speed of conver- 
gence depends on the source size. The convergence in terms of 
peak flux density and integrated flux for the sources is faster 
for compact sources than for extended ones. For sources smaller 
than 320", the convergence is reached within 1% in peak flux 
density and 2% in total flux. For the few larger sources, the con- 
vergence is reached within 2% and 4% at iteration 20, respec- 
tively. These results are very similar to those obtained with the 
Cha I dataset (see Paper I). 
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Fig. 2. 870 /im continuum emission map of Cha III obtained with LABOCA at APEX. The projection type and center are the same 
as in FigQ] The contour levels are a, 2a, 4a, and 6a with a = 34.5 mJy/21.2"-beam, i.e. 3 times the rms noise level. The flux density 
color scale is shown on the right. The field of view of LABOCA (10.7') and the angular resolution of the map (HPBW = 21.2") are 
shown on the right. The red boxes in the insert are labeled like Figs.[3k-e and show their limits overlaid on the first 870 /mi contour. 
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Fig. 3. (continued) c Field Cha3-Center. 
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Fig. 4. a 870 /itn continuum emission sum map of Cha III at scale 5. The contour levels are -a (in dotted blue), a, 2a, and 4a, with 
a = 34.5 mJy/21.2"-beam, i.e. about 3 times the rms noise level, b Smoothed map, i.e. residuals, at scale 5. The contour levels are 
—c (in dotted blue), c, 2c, 4c, and 8c, with c — 8.1 mJy/21.2"-beam, i.e. about 3 times the rms noise level in this map. The greyscales 
of both maps are different. The sum of these two maps is strictly equal to the original map (Pig. |2j. 
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Fig. 6. Distributions of physical properties obtained for the 29 starless sources found with Gaussclumps in the sum map of Cha III 
at scale 5. The mean, standard deviation, and median of the distribution are given in each panel. The asymmetric standard deviation 
defines the range containing 68% of the sample, a FWHM sizes along the major (solid line) and minor (dashed line) axes. The 
filled histogram shows the distribution of geometrical mean of major and minor sizes. The mean and median values refer to the 
filled histogram. The dotted line indicates the angular resolution (21.2"). b Same as a but for the deconvolved sizes, c Aspect ratios 
computed with the deconvolved sizes. The dotted line at 1.4 shows the threshold above which the deviation from 1 (elongation) 
can be considered as significant, d Peak H2 column density. The dotted line at 5.0 x 10 21 cirr 2 is the 5cr sensitivity limit, e Visual 
extinction derived from 2MASS. 
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Fig. 7. Distributions of masses and free-particle densities obtained for the 29 starless sources found with Gaussclumps in the sum 
map of Cha III at scale 5. The mean, standard deviation, and median of the distribution are given in each panel. The asymmetric 
standard deviation defines the range containing 68% of the sample, a Peak mass of the fitted Gaussian, b Mass within an aperture 
of diameter 50". c Total mass of the fitted Gaussian. The dotted line indicates the estimated completeness limit at 90% for Gaussian 
sources corresponding to a 6.3<x peak detection limit for the average source size, d Mass concentration, ratio of the peak mass to 
the mass within an aperture of diameter 50". e Peak density, f Mean density within an aperture of diameter 50". g Mean density 
derived from the total mass, h Density contrast, ratio of peak density to mean density. In panels a and e, the dotted line indicates 
the 5<x sensitivity limit. The upper axis of panel d, which can also be used for panel h, shows the power-law exponent p derived 
assuming that the sources have a density profile proportional to r~ p and a uniform dust temperature. Alternately, the upper axis of 
panel h, which can also be used for panel d, deals with the case where the density is uniform within a diameter Df\- dt and decreasing 
as r~ 2 outside, still with the assumption of a uniform temperature. 
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Fig. 8. a Total mass versus mean FWHM size for the 29 starless sources found with Gaussclumps in the sum map of Cha III at 
scale 5. The angular resolution (21.2") is marked by the dotted line. The solid line (M oc FWHM 2 ) is the 5cr peak sensitivity 
limit for Gaussian sources. The dashed line shows the 6.3<x peak sensitivity limit which corresponds to a completeness limit of 
90% for Gaussian sources, b Total mass versus mean deconvolved FWHM size. Sizes smaller than 25.4" were set to 25.4" before 
deconvolution (see note b of Table [3]). The solid line shows the relation M = 2ARa 2 /G that characterizes critical Bonnor-Ebert 
spheres (see Sect. 14.3.61) . The dash-dotted line shows the location of this relation when divided by 2, and the dashed lines when 
divided by 4 and 6. The source with a mass larger than the critical Bonnor-Ebert mass Mbe,f„, estimated from the ambient cloud 
pressure is shown with a filled circle, c Mass concentration versus mean FWHM size. The dashed line is the expectation for a 
circular Gaussian flux density distribution. 
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Fig. A.l. Statistical properties of 100 artificial, circular, Gaussian sources inserted into the raw time signals before data reduction (a 
to e) or directly inserted into the final unsmoothed continuum emission map of Cha I (f to j). a and f : ratio of output to input peak 
flux density as a function of input peak flux density, b and g: ratio of output to input peak flux density as a function of input size 
(FWHM). c and h: ratio of output major size to input size as a function of input size, d and i: ratio of output minor size to input size 
as a function of input size, e and j: ratio of output to input total flux as a function of input size. The size of the black symbols scales 
with the input peak flux density of the sources. In addition, crosses and plus symbols are for input peak flux densities below and 
above 150 mJy/19.2"-beam, respectively. The red circles in panels b to e and g to j show a fit to the data points using the arbitrary 
3 -parameter function y = log(a/(l + (x/fiy)). 
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Fig. A.3. Convergence of the iterative data reduction for artificial, circular, Gaussian sources with a peak flux density larger than 150 
mJy beam 1 . The first and second columns show the ratio of the output to the input peak flux densities, and its relative variations, 
respectively. The third and fourth columns show the same for the total flux. The different rows show sources of different widths, 
as specified in each panel of the first column. The size of the symbols increases with the input peak flux density. In panels n and 
p, the noisier curve corresponds to a very large source with a peak flux twice as weak as the other one that has a similar size (see 
Figs. IA. la and b). 
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